Mon. Not. R. Astron. Soc. OOP. [TlfTSl C2008') Printed 8 May 2008 (MN WI^ style file v2.2) 



Hydrodynamical simulations of Galactic fountains I: 
evolution of single fountains. 



C. Melioli,^* F. Brighenti\ A.D'Ercole^ and E.M. de Gouveia Dal Pino^ 

^ Dipartimento di Astronomia, Universita di Bologna, via Ranzani 1, 40126 Bologna, Italy 
^ INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 40126 Bologna, Raly 
^lAG, Universidade de Sao Paulo, Rua do Matdo 1226, 05508-090 Sao Paulo, Brazil 



Accepted ... Received in original form ... 



ABSTRACT 

The ejection of the gas out of the disk in late-type galaxies is related to star for- 
mation and is due mainly to Type II supernovae. In this paper we studied in detail 
the development of the Galactic fountains in order to understand their dynamical 
evolution and their influence in the redistribution of the freshly delivered metals over 
the disk. To this aim, we performed a number of 3D hydrodynamical radiative cool- 
ing simulations of the gas in the Milky Way where the whole Galaxy structure, the 
Galactic differential rotation and the supernovae explosions generated by a single OB 
association are considered. A typical fountain powered by 100 Type II supernovae may 
eject material up to ~ 2 kpc which than collapses back mostly in form of dense, cold 
clouds and filaments. The majority of the gas lifted up by the fountains falls back 
on the disk remaining within a radial distance AR = 0.5 kpc from the place where 
the fountain originated. This localized circulation of disk gas does not influence the 
radial chemical gradients on large scale, as required by the chemical models of the 
Milky Way which reproduce the metallicity distribution without invoking large fluxes 
of metals. Simulations of multiple fountains fuelled by Type II supernovae of different 
OB associations will be presented in a companion paper. 
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1 INTRODUCTION 



Deep H I observations and a detailed mode lling of the H 
I in edse-on spiral galaxi es as NGC 891 llSwaters et al. 
19971 : lOosterloo et al.il2007l ). UGC 7321 l|Matthews fc Wood 



20031) or part i ally i n clined spiral galaxies as NGC 2403 



Schaap et all l2000l: iFraternali etafl |2002| ). NGC 4559 



Barbieri et al.l 20051 ) have shown the presence of thick 



gaseous layers extending to several kiloparsec from the plane 
of the disk, with a mass of 1/10 of the total H I. Optical ob- 
servations have also rev ealed extraplanar emission from dif- 
fuse ionised gas (DIG) (iFerguson et al.|[l996l : iHoopes et al.l 
1 19991 ) which, as the H I, extends to a large extent above 
the plane. The extraplanar gas (EPG) has peculiar kine- 
matics that has been studied both for the neutral and the 
ionised component. The modelling of the H I layer has shown 
that the EPG rotates in the same direction as the disk, 
and has a mean rotation velocity about 20-50 km s~^ lower 
than that of the gas on the plane (e.g. ISwaters et al.|[l997l : 
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IFraternali et al. 2004). S imilar results hold for the ionised 
component (jFraternali et al. 2004). 

In general, the extraplanar gas appears to be regularly 
distributed over the whole galaxy , and its amount is related 
to th e rate of star formation (S chaap et al.|[2000l : iDettmail 
|2005| . e.g.). This suggests that a significant fraction of the 
EPG has an internal origin, likely being disk ISM ejected in 
the halo by (clustered) supernova explosions. 

The presence of gi ant holes in the H I d isk of many 
galaxies as NGC 6946 (iBoomsma et al.ll2005l). NGC 2 403 
llFraternah et al.ll2002l), M3 3 dDeul fc den Harto3ll99d) , M 



101 ("Kamphu is et al.l ll99l'). M31 (Brinks & Baiaia '198' 
and in the Milky Way by (Hcilcs 1984) and (Koo ct al. 199: 
also point to an internal origin of the EPG. These holes 
can be hardly created b y the competing scenario of infalling 
clouds (|Comer6nll2003 ). Instead, they can easily accounted 
for by associations of O and B stars exploding in a corre- 
lated fashion and giving rise to superbubbles. In some cases 
vertical motions are associated with these h oles, as expected 
once the superbubble brea ks out of the disk (|Mac Low et al.l 
ll989l : lKoo fc McKed [l993 ) . reinforcing the hypothesis that 
they are linked to galactic fountains. 
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The above arguments point at a substantial gas circu- 
lation between disk and halo. The ejection of the gas out 
of the disk is thought to be related to star formation in 
the disk and to be due mainly to Type II supernovae (SNe 
II). This process is usually described as galactic fountain 
ijShapiro fc Field! 1 19761 : iBregmanI [l980l : lKahnlll98ll l. In the 
simple version, the scenario assumes that gas heated by SNe 
explosions rises above the plane of the disk but remains 
bound to the galaxy. Then, the material will recirculate, 
falling backward to the galactic plane as cold gas, presum- 
ably concentrated in clouds formed through thermal insta- 
bilities. According to this model, the galactic fountain gas 
moves to larger galactocentric distances as it rises above the 
galactic plane, to conserve its angular momentum. This ef- 
fect would account for the observed lower rotation velocity 
of the extraplanar gas. When the clouds form via thermal in- 
stability, they fall back toward the galactic plane and their 
radial distance from the centre decreases. From a simple 
analysis it can be shown that after about 10® yr of its forma- 
tion a cloud returns to the Galactic plane at a Galactocentric 
radius not very different from its initial value (Spitzcr 1990). 
In the framework of a one-dimensional model Kalm & Brett 
l|l993h showed analytically that, at the Sun distance from 
the Galactic centre, a cloud collide with the Galactic mid- 
plane at a distance of 850 pc from its initial position. In 
these admittedly idealized scenarios, one expects the early 
outflow phase to be mostly traced by ionised gas and the 
late inflow phase by neutral gas. 

Despite the arguments reckoned above for an internal 
origin of most extraplanar gas, a scenario in which the halo 
cold and warm gas is a ccreted from the i ntergalactic medium 
is not ruled out (e.g. iBlitz et al.lll999l ). In fact, the above 
analysis could be reversed to conclude that a higher star for- 
mation rate might be the final result of accretion rather than 
the cause at the origin of the EPG. Indeed, extended infall 
of external gas on galactic disks is needed in current galac- 
tic evolution models as a way to maintain star formation 
rates similar to the observed present ones in normal spiral s 
(e.g. iTakeuchi fc Hirashita|[200ol : ISemelin fc Combed |2002| ). 
as well as to explain the obse rved patterns in the chemi- 
cal ev olution of the Milky Way (|Chiappinill2002l : iGeiss et al.l 
|2002| ). The observational counterpart of this gas accretion 
might be the high velocity clouds (HV Cs) observed in our 
and other galaxies (|Wakkeill200ll . l2004l ). However, the abil- 
ity of these clouds to provide the needed amount of gas at 
low metallicity, as well as their capability to trigger star 
formation is still unclear. Evidence for star formation in- 
duced by HVC impacts on galactic disks is unclear, while 
numerical simulations seem to indicate th at infall-induced 
star formation is substantially prevented (|Comer6nl [iooi . 
and references therein). 

The galactic fountain scenario can also explain the dif- 
fuse soft (< 1 keV) X-ray emissio n that has been detected i n 
some spiral galaxies li ke MlO l ijSnowden fc Pietschlll995h . 
M51 jEhle et al.lll995l), M 33 (iLong et al.lll996l). NGC 89 1 
l|Bregman fc Houcklll997l). NGC 4631 jWang etHI I2OOII) 
NGC 2403 l|Fraternah et al. 1'2002'). NGC 3556 ( |Wang et al.l 
l200j). and NG C 2613 jlTet al. 2006) (for a short review 



2841 the diffuse X-ray emission arises primarily from the in- 
ner galactic disk, and its morphology is the result of various 
hot plumes sticking out of the disk and into the halo of the 
galaxy (|Warig.2007i '). These X-ray emitting plumes extend 
vertically up to a few kiloparsecs and are most likely blown 
out hot gas heated in massive star forming regions in the 
disk. This conclusion is reinforced by the observation that 
the X-ray luminosity correlates with the star formation rate 
and that the hot gas seems oxygen rich, a signature of SNII 
enrichment. 

From all the above, it turns out the fountains may be 
crucial in characterizing X-ray, as well as H I and Ha emis- 
sions in normal spiral galaxies. Moreover galactic fountains 
could also play a role in the chemical evolution of spiral 
galaxies, especially in shaping the abundance gradients in 
the disk. In principle the fountains might be able to move 
SN II ejecta relatively far from the place where they have 
been produced, affecting the gradients of a-elements such as 
Mg and O. 

The richness of phenomena involved in the galactic 
fountains can be hardly taken into account by the simple sce- 
nario outlin ed above dShapiro fc Fi eld 1976; Bregman 198^; 
lKahnlll98ll ). Nearly twenty years after these three seminal 
papers, however, the advent of powerful computers made 
possible to simulate the fountains rather accurately. 

In fact, extensive work on superbubble/chimneys and 
fountains formation and on their effects upon the inter- 
stellar medium in the disk and the halo has been per- 
formed over the last decades (e.g. Tomisaka fc Ikeuch 
'1986; Norman & Ikcuchi 1989; Mac Low & M cCray 19&I. 
Ma c Low et all 1 19891 : iHouck fc: Bregmanlll990i : iKorpi et al] 
I1999i r most of which has shown 2D and 3D simula- 
tions considering the evolution of bubbles and super- 
bubbles both in unmagnctizcd and magnetized media. 

In particular, iTenorio-Ta glc fc Bodcnhcimcr ( 19 8§i) and 
' Norman fc Ikeuchil (Il989l) (see also iTomisaka fc Ikeuchil 



m 



mn3|2003)). This hot gas does not extend significantly 



more than 10 kpc away from the galactic disks, and appears 
to have substantial substructures which are best appreci- 
ated in disk galaxies that are moderately inclined. In NGC 



1 19861 : Mac Low et al.l 19891 ). were the first to explore the 
role that superbubles and chimneys physically connected to 
underlying OB associations in the disk play o n the foun- 
tains evolution in galaxies. iHouck fc BregmanI ( |1990 ) bor- 
rowig concepts from earlier works, investigated the physical 
conditions for the formation of subsonic and transonic foun- 
tains as a function of the rad i ative cooli ng time. 

More recently, Ide Avillej (|200G| ) and|de Avillez fc Berrvl 

(l200ll ) performed a number of 3D hydrodynamical simula- 
tions of the gas in the Milky Way in order to account for 
the collective effects of supernovae on the structure of the 
interstellar medium (ISM). In order to obtain a high spatial 
resolution, only a small region of the Galaxy has been con- 
sidered by the authors; the simulated volume has an area of 
1 kpc^ and a vertical extension varying between 4 and 10 
kpc. The gas in this volume is initially assumed to be plane 
parallel stratified, and any dependence of the variables on 
the Galactocentric distance is neglected. The general struc- 
ture of the ISM is reproduced, with a thin H I disk overlaid 
by thick H I and H II gas disk. This structure is essentially 
given by SNe II isolated or assembled in small clusters. The 
gas powered by SNe II associated in larger clusters reaches 
larger heights, where it cools forming clouds which eventu- 
ally descend. Intermediate velocity clouds (IVCs) form pref- 
erentially in the range 0.8 < z <2 kpc whereas HVCs (10% 
of total number of clouds) form at larger heights. 
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Three-dimensiona l simulations have been reported also 
bv iKorpi et all 

who include the galactic magnetic 
field in the disk and the Galactic differential rotation, but 
the vertical extension of the simulated volume was ± 1 kpc, 
covering an area of 500 by 500 pc^ with a resolution of 8 
pc. The lack of an extended grid made the development of 
a disk-halo-disk cycle impossible, because part of the mass 
is lost through the grid boundaries. MHD sim ulations were 
carried out also bv lde Avillez fc Breitschwerdd ((2005.), with- 
out including differential rotation, where again in order to 
reach very high resolution (of 0.6 pc) they have considered 
only a small volume of the Galaxy. These MHD studies have 
revealed that the gas transport into the halo is not prevented 
by the parallel magnetic field of the Galaxy, but only delayed 
by f ew tens of Myr when c ompa red to pure HD simulations. 

iFraternali fc BinnevI (|2006l ) followed a different strat- 
egy. They considered the whole Galaxy, and calculated or- 
bits of gas clouds depicted as bullets that are ejected (by 
SNe II) from the disk and then move ballistically up into 
the halo and back down to the disk. Hydrodynamical effects 
are essentially neglected in these models. 

In this paper we would like to study the motion of 
the gas on large scales in order to understand the ob- 
served kinematic of the EPG and the possible influence in 
shaping the metallicity gradient in the disk. To this aim, 
we need different kind of simulations with respect to the 
aforementioned papers |Korpi et al. I ll999l : lde Aville j l2000l: 
Ide Avillez fc Berrvf I2OO1I ). In fact, although iKorpi et all 



119981) take into account the Galactic rotation, the simulated 
volume is to small to follow the gas circulation over large 
distances. T he volumes co nsidered bv lde Aville j l|2000l ) and 
Ide Avillez fc Berrvl l|200ll ) are larger, but in this case the 
rotation is neglected and any radial depe n dence is absent. 
On the other hand, IFraternali fc BinnevI l|2006l ) take into 
account the whole Galaxy, but their approach is purely bal- 
listic. Our approach is to consider the whole Galaxy and run 
3D hydrodynamical simulations of the fountains. While the 
numerical grid used here includes the whole rotating Galaxy, 
our adaptive mesh scheme allow us to follow in detail only 
a limited volume of it, the remaining galactic volume being 
mapped at a lower resolution. Although a detailed descrip- 
tion of the different phases of the ISM is hampered by this 
limited resolution, the large scale dynamic and thermal evo- 
lution of the fountains can be satisfactorily followed. Here 
we present simulations of single fountains, i.e. generated by 
a single OB association. We shall describe simulations of 
multiple fountains in a companion paper. 



2 THE MODEL 
2.1 Galaxy model 

The ISM in our model is made up of three components, 
namely molecular (H2), neutral (H l) and ionized (H II) hy- 
drogen. Following I Wolfire et al.l l|2003l ), each density compo- 
nent in the disk is assumed to be of the form: 



p = — exp 

2zd 



R_ 

Rd 



Hi 

Zd 



(1) 



where R is the cylindrical radius, Zd is the vertical scale 
height and Rd is the radial scale length of the disk. The pa- 
rameter i?m allows for the depression in the gas density ob- 



served in the inner several kiloparsecs of the Galaxy. Ed rep- 
resents the superficial density of the different components. 
The values of the parameters are summarized in Table [T] 

The ISM is initially set in rotational equilibrium in the 
Galactic gravitational potential given by the summation of 
dark matter halo, bulge and disk contributions. 

The dark matter halo gravitational potential is as- 
sum ed to follow th e Navarro, Frenk and White profile 
( Navarro et al.ll 19961 ) 



^ r \ A 2 ln(l -I- x) 

9dni{r) = -4-7rGrdni,0/9dm,0 



(2) 



where pdm,o is a reference density, rdm.o is a scale radius, 
X = r/rdm,o and r is the spherical radius. The halo is trun- 
cated at a radius rdm,t beyond which its potential follows 
the 1/r profile. See Table 2 for the numerical values of these 

parameters. 

The bulge gravitational potential is given by (|Hernauistl 

Il990f ): 



rh,o + r ' 



(3) 



where rb,o is a scale radius and Mb is the bulge mass. 

Finally, the gravitational potential of the disk is as- 
sumed to be generated by a stellar distribution following 
a flattened King profile 

P*,o 



[1 + {R/R.,cf + {z/z,,,f] 



213/2 • 



(4) 



Here p*_o is the central density of the stars, and _R*_c and 
z*,c are the core radii, whose ratio is 5 = 2;*,c/J?*,c = 
0.03 In order to avoid an unbounded growth of the stellar 
m ass with radius the st ellar profile is truncated wherever 

yJ{R/R-,,cY ^ (2/z*,c)2 > R*,t/R.,c = Z*,t/2*,c, with R^,t 

and 2*,t tidal lengths whose ratio is z*,t/i?*,t ~ 5. The stellar 
potential is computed numerical ly fol lowing the method 
described in Brighonti & Mathcw^ (|l996l ). 

Table Ogives the values of all the parameters concerning 
the Galaxy model. With such values we obtain Mdm = 1-3 x 
10^^ Mq, Mb = 3.5 X 10^" M© and PL = 9.6 x 10^" Mq 
for the dark matter halo, the bulge and the stellar disk, 
respectively. 

In order to set the ISM given by Eq.[T]in a rotating con- 
figuration in equilibrium in the general potential well of the 
Galaxy, we proceed as follows. The pressure at any point is 
found integrating the z-component of the hydrostatic equi- 
librium equation for any value of the disk radius R. The 
integration starts at the outermost values of z (where we 
can assume the pressure P = 0) and proceeds inward to 
reach the galactic plane z — 0. 

We then obtain the rotation velocity from 



vliR, z) = vl{R, z) - — 
P 



dP(R,z) 



dR 



(5) 



where Vc{R,z) — ^ R{d${R, z)/dR)z=cost is the circular 
velocity, and $ is the total potential of the Galaxy: $ = 
'^dm + $b + The values v^{R,0) obtained in this way 
mimic quite well the rotation curve of the Galactic disk (c.f. 

Fig. ID). 

As shown in Fig. O that contains the temperature, 
pressure and density distributions of the Galaxy, the gas 
has a temperature always greater than 10'* K. This comes 
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from two conditions. On one hand, from the necessity to 
hold in equil ibrium the disk gas d istribution against gravity, 
as given by IWolfire et all (|2003h . On the other hand, and 
more important, it is due to the assumption that H I has a 
smooth (rather than a clumpy) distribution, so that a ther- 
mal, rather than a turbulent, pressure is required to balance 
gravity. 

In most models a hot (Th = 7 x 10® K) isothermal gas 
halo is added which is in equilibrium with the Galactic po- 
tential well: 



/9h = Po,h exp[- z) ~ e^^{R, 0) - (1 



-e^)$(0,0))/c^] 



(6) 



Here po = 2.17 x 10~ g cm~ is the halo central density, 
and c is its isothermal sound speed. In some models the 
hot halo is allowed to rotate with velocity u^^h ~ evc{R,0), 
where e is an arbitrary parameter with values in the range 
e ^ 1. In the models where the halo is present, the disk 
density is replaced by the hot halo density wherever the halo 
pressure is larger than the pressure of the disk. When this 
halo is included, the initial pressure values which are needed 
to integrate the z-component of the hydrostatic equilibrium 
equation (see above) are given by the halo pressure at the 
z-boundary. 



Figure 1. Galactic rotation curves for different values of z. 
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Figure 2. Edge-on logarithmic distribution of temperature (top 
panel), pressure (middle panel) and density (bottom panel) of 
the unperturbed ISM in our Galaxy model. Distances are given 
in kpc, temperatures in K, pressure and density in c.g.s. 



Table 1. ISM parameters. 
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H II 


1.39 





30.0 
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2.2 Supernovae explosions 

It is known that th e great majority of th e massive stars 
form in clusters (e.g. I Ward-ThompsonI |2002| . and references 
therein). In the following, we define "single fountain" as a 
fountain powered by the SNe II present in a localized sin- 
gle cluster (for the sake of simplicity, with SNe II we mean 
here all the explosions due to core collapsed stars, i.e. gen- 
uine SNe II, S Ne lb and SNe Ic). At the observed rate of 
1.2 X 10"'^ yr l|Cappellaro et ai1ll997l ). several single foun- 
tains may occur sufflciently close to each other in time and 
space to interact mutually. This interaction obviously mod- 
ifies the general dynamics of the SN II ejecta of a single 
fountain. In order to understand the effect of different fac- 
tors (such as the Galactocentric distance, the presence of a 
gaseous halo, the rotation velocity of this halo), and how 
these effects are influenced by the interaction among several 
fountains, we have computed two classes of models differing 
only in the assumed number of stellar associations, i.e., sin- 
gle fountain models (SFMs) and multiple fountains models 
(MFMs). Here we focus on the first class of models, while 
the latter will be discussed in a forthcoming paper. 

2.2.1 Single association 

In the case of a SFM, we assume that the OB association is 
located on the Galactic plane {z = 0), and is totally housed 
inside a single mesh of the highest resolution level of our 
adaptive grid (cf. sect. [3]). The SNe II do not explode all 
at the same time because of the different masses of their 
progenitors. Thus the energy and the mass delivered by all 
the explosions are released over a time interval At=30 Myr 
(the lifetime of an 8 Mq star, which is the least massive 
SN II ancestor). We inject the explosion energy in the form 
of thermal energy at a rate E = TZEq, where Eo = 10^^ 
erg is the energy of a single explosion, TZ = A^sN/Af is the 
SN II rate, and A^sn is the total number of SN II explo- 
sions occurring in the association. The mass of the SN II 
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ejecta is injected at a rate M — TZMcj , where Mej = 16 Mq 
is the mean mass expelled by a single SN II explosion. As 
a result, we do not consider a sequel of instantaneous ex- 
plosions but, instead, a wind inflating a bubble with a me- 
chanical luminosit y Lw = E. This treatment is supported 
by the findings of iMac Low fc McCravl () 19881 ') who showed 
that discrete explosions can be approximated as continuous 
events when more than ~ 7 SNe co me one after the other . 
About the temporal wind evolution, iLeitherer et al.l (|l999l ) 
showed that SNe belonging to an istantaneous burst release 
mass and energy at rates oscillating around a mean value; 
the general properties of the fountain we are looking for de- 
pend on such mean values rather than on the oscillations. 
We thus assumed constant values of M and Lw, as described 
above. 

As the bubble expands, the swept up ISM is compressed 
in a dense radiative shell edging the bubble itself. The inte- 
rior of the bubble, instead, is filled by the hot, low density 
ejecta (c.f. lWeaver et al.lll977l ). The bubble expands faster in 
the direction perpendicular to the plane, because the den- 
sity decreases more steeply in that direction. If the den- 
sity gradient is strong enough, the shell can accelerate. The 
acceleration of the shell causes it to be disrupted by the 
Rayleigh- Tailor instability, and the hot gas inside the bub- 
ble subsequently leaks out and drives a new shock into the 
ambient medium, giving rise to the galactic fountain. The 
break out is expected to occur wh en Lw > 2.5I/b, where the 
critical luminosity is given by Cc.f. lKoo fc McKeelll992h . 



Lb = 2.4 X 10^^Po,4i?cff,2Co,6 erg s'V 



(7) 



Here Po,4 = Po/{W*k cm~"^ K), where Po is the midplane 
pressure, and k is the Boltzman's constant. co,6 = co/(10^ 
cm s~^), where cq represents the sound speed in the mid- 
plane; actually, as the temperature varies somewhat within 
the range < z < H^s, we adopt the mean value of c in 
Eq. [7] Finally, the effective scale heigh of the ISM at the 
Galactocentric distance R is defined as 



p{R,0) 



p{R, z)dz 



(8) 



All the SFMs we run have the same number A^sn = 100 su- 
pernova explosions. This number h as been chosen because 
typic al for Galactic OB associations (|Higdon fc Lingenfelteij 
I2OO5I . see also Paper II) . Furthermore, such a number assures 
the realization of the break out at any location of the foun- 
tain on the disk (cf. sect. 15.1) ). 

It is worthwhile to note that our numerical grid is at 
rest in space, while the Galaxy rotates. Thus, the OB asso- 
ciation moves around describing a circle on the z = Q plane. 
In order to properly locate the energy and mass source in 
the grid, we compute separately the trajectory of the asso- 
ciation, and at each time we put an appropriate amount of 
thermal energy and mass in the smallest cell of our multi- 
level grid (see below) transited by the association at that 
time. 



3 THE NUMERICAL METHOD 

In our simulations we have used a modified version of the 
adaptive mesh refinement YGUAZU code. This code inte- 
grates the 3D inviscid gasdynamic equations (expressed in 



Cartesian coordina t es) w ith the flux vector splitting algo- 
rithm of Ivan Leer! l|l982l ). Radiative losses are computed 
through a non-equilibrium ionization (NEI) calculation in- 
cluding a set of continuity equations for atomic/ionic or 
chemica l species; in the pr esent simulations solar abun- 
dances (|Asplund et al.l |2005| ) are assumed. The details of 
the gasdynam ic and the adaptive grid a lgorithms have bee n 
presented by iRaga et all (|20oO) and iRaga et all (|2002h . 
and tests of the code, as well as, s everal physical applica- 
tions can be found in i Raga et al.1 ||2002|). iMasciadri et al.l 
(|2002h . iGonzalez et al.1 (|2004| ) "and iMelioh et al.1 \2Qm . 
When compared with other methods, this particular scheme 
has a reasonable performance (as far as diffusion is con- 
cerned) when one consid e rs, e.g., ID test problems (se e, e.g., 
Ivan Albada et al.1 1| 19821 ): I Woodward fc Coleilal l|l984l )'). Ar- 
tiflcial viscosity is set to zero in all our simulations. As dis- 
cussed in the previous Section, at any time the mass and 
energy source terms due to the SNe are absent everywhere 
but in the cell hosting the OB association. 

The 3D binary, hierarchical computational grid is struc- 
tured with a base grid, and with a number of nested grids 
whose resolution doubles going from one level to the next 
one. In our most refined simulations, we have adopted six 
grid levels (with the linear size of the meshes of the most 
refined grid being 13 pc) covering a box with physical di- 
mensions of 26.4 X 26.4 x 6.6 kpc. The grid being adap- 
tive, in principle the maximum resolution should be set at 
the boundary between the hot halo and the disk, where the 
density gradient (along the z direction) is very steep. This, 
however, would represent an useless waste of cpu time be- 
cause we are interested in describing in particular the foun- 
tain which evolves in the selected area of the disk. We thus 
enforce the maximum grid resolution only in this area and 
in the volume above it, even if the density gradients in the 
halo are mild. In fact, we want to follow the circulation and 
the thermal history (i.e. the degree of radiative cooling) of 
the metals expelled by the SNe II, as carefully as possible. 
For this reason we have added three different tracers pas- 
sively advected by the code describing the disk gas, the halo 
gas and the SN II ejecta. In this way, we can distinguish 
the evolutive history of each gas component and understand 
their degree of mixing, as the fountain evolves. 

Given the symmetry of the problem, our computational 
volume encompasses only the "upper" half of the Galaxy. We 
thus have outflow boundary conditions everywhere but at 
"the bottom" (i.e. the z = plane) where reflecting bound- 
aries are enforced. 



4 THE REFERENCE MODEL RM 

In this section we discuss extensively the properties of the 
RM, which we take as a reference model. In this model the 
gaseous halo is present but does not rotate, and the OB as- 
sociation is located at the galactocentric distance R = 8.5 
kpc. At this radius the transition between the disk and the 
hot halo occurs a,t z — 800 pc. The critical luminosity is 
Lb = 1.5 X 10^^ erg s~^, well below the mechanical luminos- 
ity Lw = 10'^* erg s~^ provided by the 100 SNe II powering 
the fountain. This model has been run with a maximum 
spatial resolution of 13 pc. 

In Fig. |3] we show the face-on column density of the SN 
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Table 2. Galactic parameters. 
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Figure 4. Each row represents the face-on view distribution of the z column density of the SN II ejecta up to 2=400 pc (left panel), 
2=800 pc (central panel), and z =6.6 kpc (right panel) for RM; the upper, middle and lower rows refer to t =30 Myr, t=114 Myr and 
t=160 Myr, respectively. The circles visible in the panels have the same meaning than in Fig.|3] The logarithmic column density scale is 
given in cm^-^ and the x — y scale is in kpc. 



II ejecta after t = 160 Myr. During its activity the fountain 
digs a hole on the disk, and throws SN II ejecta and ISM 
vertically up to |z| = 2 kpc above the galactic plane. Once 
the stellar explosions cease, the hole collapses in ~ 2 x 10^ 
yr and the ejecta trapped at its edges (nearly half of the 



total, see below) mixes with the local ISM. Owing to the 
differential rotation of the Galactic disk, the ejecta does not 
remain confined in one point, but is stretched, giving rise to 
the bean- like structure seen in Fig. |3l The low density "tail" 
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Figure 3. Face-on view of the z column density distribution of 
tiie SN II ejecta after t = 160 Myr for RM. The SNe II explode at 
a radius R = 8.5 kpc, halfway between the radii R = 8 kpc and 
R = 9 kpc of the two circles drawn in figure. These two circles 
are drawn to guide the eye, and their distance represents approx- 
imately the maximum extension of the hole in the ISM carved by 
the fountain on the disk during its activity. At the time to which 
the figure refers to the hole has collapsed and disappeared. The 
logarithmic column density scale is given in cm"'^ and the x 
scale is in kpc. 



y 



with a banana-shape is given by the ejecta pushed at high 
altitude which then slowly falls back lasting above the disk. 

The two concentric circles in Fig. [3] with radii R = 8 
kpc and R = 9 kpc, have been drawn to guide the eye. 
The SNe II explode between these two circles, and their dis- 
tance corresponds to the diameter of the hole produced by 
the fountain. Comparing the shape of the tail traced by the 
ejecta with the circles, we note that the gas of the fountain 
tends to move inward during its trajectory, rather than out- 
ward, as one could expect, at least initially (cf. section [l]). 
Actually, as the gas moves upward and interacts with the 
halo, it transfers to it part of its angular momentum; the 
centrifugal force decreases in pace with the circular veloc- 
ity and the centripetal component of the gravity prevails, 
pushing the gas toward the Galactic centre. 

To better illustrate the fountain evolution, we show in 
Fig. |4] a sort of tomography of the fountain itself. Each 
single row in this figure represents the column density of 
the ejecta up to 2=400 pc (left panel), z=800 pc (central 
panel), and z >6.6 kpc (right panel); the upper, middle and 
lower rows refer to t =30 Myr, t=114 Myr and t=160 Myr, 
respectively. From the first row, we see that, initially the 
gas pushed by the fountain rises quite vertically forming a 
column whose base has a size comparable to the distance 
between the two circles. In the remaining panels the pro- 
gressive "stretching" of the ejecta is evident, as well as the 
formation of the "cometary" tail at high values of z and its 
tendency to bend toward the Galactic centre. 

A further insight on the evolution of the galactic foun- 
tain is obtained from Fig. [5] The three panels illustrate the 
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Figure 5. The upper panel shows the evolution of angular mo- 
mentum of the SN II ejecta for the RM. The evolution of the 
amount of the ejecta above different heighs is illustrated in the 
middle panel. Finally, the lower panel displays the temporal be- 
haviour of the amount of the ejecta located at i? < 8 kpc (dashed 
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evolution of several quantities starting from t=30 Myr, the 
time at which the SNe II stop to explode. The upper panel 
shows how quickly the ejecta loses its angular momentum 
because of the interaction with the gaseous halo. After 80 
Myr nearly 10% of the angular momentum of the fountain 
has been transferred to the hot halo; later on, no further 
transfer occurs because nearly 75% fraction of the ejecta is 
located below the disk-halo transition and rotates toghether 
with the ISM Q. This is illustrated by the middle panel of 
Fig. [5] showing the amount of the ejecta mass fraction lo- 
cated above different heights. It is interesting to note that, 
apart for some fluctuations, for t >80 Myr, the long-term 
evolution of these quantities becomes very slow. This is due 
to the fact that most of the ejecta situated above the plane 
is rather diluted, and tends to float together with the extra- 
planar ISM. Finally, in the lower panel of Fig. [5] we quan- 
tify the tendency of the ejecta to move radially by plotting 
the fraction of the total mass of ejecta located at R >9 
kpc and R <8 kpc (i.e. the radii of the two circles drawn 
in Figs. |3] and |4)|; the amount of mass located within the 
region 8 < 7? < 9 kpc is not taken into account. In the be- 
ginning the ejecta starts to follow the expected tendency to 
move outward but, as the loss of angular momentum pro- 
ceeds, the fraction of gas moving inward increases and after 
60 Myr overrides that directed outward. The amount of the 
outer mass remains always very low and does not evolve very 
much; on the contrary, in the inner region the mass of the 
ejecta keeps increasing during almost all the period of the 
simulation. 

In Fig. [5] are also shown results obtained with a coarser 
spatial resolution (of the finest grid) of 26 pc. Actually, in 
order to test the convergence of our simulations, we have run 
simulations with different higher refinement level resolution: 
52, 26 and 13 pc. The gross features of the hydrodynamical 
evolution of the fountain over the whole gas cycle do not 
change significantly and are captured even in the coarser 
case. We notice that in the coarser case there is a slightly 
larger deposition of material at higher latitudes (10% larger) 
than in the higher resolution case due to smaller radiative 
cooling. As a consequence, the coarser case has also a slightly 
larger loss of angular momentum, however the overall behav- 
ior of the fountain evolution is found to be very similar in 
both cases and does not affect the main conclusions of the 
present work. 

The gas lifted up by the fountain has a mass 2.5 x 10^ 
M0 , almost all (92%) condensed into dense filaments cooled 
to r = 10'' K, the minimum allowed temperature (see Fig. 
[6] and Fig. [7]). The clouds form via thermal instabilities at 
z ~ 2 kpc once the ascending gas has inverted its motion 
and moves back toward the disk. In fact, the clouds have all 
negative velocities in the range 50-100 km s"\ We point 
out that the clouds in our simulations cannot form by the 
hot wind gas powered by the SNe II. In fact, the sound 
crossing time of this gas is shorter tha n the cooling time, 
and cold condensations are prevented l|Houck fc BregmanI 




R (Kpc) 

Figure 6. Vertical logarithmic distribution of the gas at t = 50 
Myr on a plane orthogonal to the Galactic plane and passing 
through the place where the SN II association is located. Log- 
arithmic isodensity contours (g cm~^) of the ejecta are overim- 
posed, highlighting the fountain pattern and the cloud formation: 
the ten contours cover the range between -28.7 and -26.7, with 
constant intervals of 0.2. As the plane in figure does not encom- 
pass the Galactic axis, the gas disk shows an apparent outward 
radial velocity given by the projection of its rotational velocity. 
This projection corresponds to 75 km s~^ 




7.4 7.8 8.3 8.7 9.7 9,6 
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Figure 7. Vertical logarithmic distribution of the gas tempera- 
ture at t = 50 Myr. Logarithmic isopressure contours (dyn cm~'^) 
of the pressure are overimposed: the fifteen contours cover the 
range between -13. and -10.8, with constant intervals of 0.146. 



^ The transfer of angular momentum is evaluated comparing the 
angular momentum of the ejecta within the grid with the origi- 
nal angular momentum held by the ejecta at the moment of its 
injection within the grid. 
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Il990l ). Instead, the clouds originate by the compression and 
cooHng of the disk gas Ufted by the fountain. 

The chemical composition of these clouds is practically 
unaffected by the SN II ejecta. The fountain is powered by 
100 SNe II, half of them exploding in the half space mapped 
by the gri d; as each supernova d elivers on average 3 M© of 
metals (cf. iMarcolini et al.ll2007h , a total mass of 150 Mq of 
heavy elements is ejected by the fountain within the compu- 
tational domain. At f = 50 Myr, when the lifted ISM reaches 
the maximum height, 50% of these metals results trapped 
within the disk, while ~ 30% remains diffuse floating over 
the disk as hot, diffuse gas and only ~ 20% is locked in the 
clouds. As a result, the metallicity increment in the clouds 
due to the freshly delivered metals corresponds to ~ 0.01 in 
solar units, negligible compared to the ~ solar abundance of 
the ISM. As we will show in a forthcoming paper, this result 
also holds in models with multiple SN II associations. In 
conclusion, almost all the gas lifted up by the fountain con- 
dense into clouds without being chemically affected. After 
150 Myr 45% of the fresh metals stays "on the disk" (below 
z = 800 pc) within a radial distance of AR = 0.5 kpc from 
the OB association. A further fraction of 35% is found on 
the disk within the range 9.5< R <7 kpc. The remaining 
20% of metals is still over the disk, half of it at _R > 8.5 kpc 
and haff at i? < 8.5 kpc. (cf. Fig. Hand Fig. i|. 



5 CHANGING THE PARAMETERS 

In this section, we illustrate the differences among the results 
obtained with our reference model and those taking different 
parameters, the aim being to understand the role played by 
these parameters in the evolution of a galactic fountain. In 
order to obtain a relatively large number of models, we have 
run them adopting a spatial resolution for the finest grid 
(26 pc) which was worse than that adopted for the reference 
model (13 pc). This latter model was also run with this 
coarser grid in order to make homogeneous comparisons. 

Although the results depend on the adopted spatial res- 
olution, several general characteristics of RM are not greatly 
affected by the minimum mesh-size adopted, as shown in 
the top and bottom panels of Fig. [5] In the middle panel, 
instead, we note that the peaks present at t = 40 Myr are 
higher, while the successive "plateau" are lower, the differ- 
ence being grater for the dashed lines. The higher peaks are 
due to the smaller volume from which the SN II energy is 
delivered: in the finer grid the same energy is in a smaller 
volume and the resulting pressure pushing the gas is larger. 
The lower plateau derives by the improved resolution with 
which the clouds are followed by the finer grid: in this case, 
in fact, clouds are smaller and denser, and can fall faster to- 
ward the disk. These differences, however, will not influence 
the relative comparison below among models. 



5.1 Varing the Galactocentric distance 

In this section, we explore the influence of the galactocentric 
distance on the fountain dynamics. In particular, we run two 
models similar to RM, but located at i? = 4.5 kpc (the inner 
fountain) and R = 10.5 kpc (the outer fountain). At these 
distances, the corresponding critical luminosities (cf. sect. 
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Figure 9. Variation of the critical luminosity with the Galac- 
tocentric distance, normalized to -Lcr,0 = 1-5 X 10^^ erg s~^, 
the value of the critical luminosity at the solar position R = 8.4 
kpc. The line represents an useful analytical approximation in the 
range of interest. 

\2:2A} are Lb = 4 x 10^'' erg s"^ and Lb = 10^^ erg s"^ re- 
spectively (the general dependence of the critical luminosity 
on the Galactocentric distance is given in Fig. |9]). 

Figure |8] illustrates the cavities carved into the ISM by 
these two fountains (upper panels) , as well as the ejecta dis- 
tribution, at t = 30 Myr. Both "chimneys" show a tendency 
to bend toward larger R because of the Galactic rotation. 
While the outer fountain experiences a well deflned break 
out, the inner fountain barely crosses the disk-halo transi- 
tion because iifcff is larger at R — 4.5 kpc. The differences 
of the effective scale heights, together with the differences of 
Lb, are also respons ible for the different shapes of the two 
fountains. Following IKoo fc McKe3 l|l992l l. it is 

i/eff V 2Tr J \mUJ ' ^ ' 

where m is the Mach number of the wind relative to the am- 
bie nt medium, and Reg is the radius of the reflected shock 
(cf IWeaver et al] ()l977h ): the contact discontinuity separat- 
ing the shocked wind from the ambient medium is at some- 
what larger radius. If Rcq/HcS is greater than unity, than 
it is possible for the wind to escape freely along the vertical 
axis without being well coUimated. If instead R^q < HcB a 
coUimated jet may form. This is actually the case for the 
inner fountain where, for m = 10, Reci/HeB = 0.85, while 
Rcq/Hcs = 1.8 holds for the outer fountain. 

More insight on the two fountains is given by Fig. 1101 
The upper panel shows the evolution of the angular mo- 
mentum of the SN II ejecta of the two fountains. The loss of 
angular momentum clearly increases with increasing values 
of R. To understand this behaviour we must consider the 
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Figure 8. Edge-on view of two single galactic fountain models. The upper panels illustrate the cavities carved into the ISM by the 
fountains at i? = 4.5 kpc (left) and _R = 10.5 kpc (right) after t = 30 Myr. The lower panels show the ejecta distributions in the two 
cases. Distances are in kpc, and the c.g.s. density and ejecta scales are logarithmic. 



diflterent ambient medium in which the fountains evolve and 
the different circulation of their ejecta. This is summarized 
in the second panel of Fig. 1101 Up to 60 Myr all the foun- 
tains push the gas essentially to the same heights (cf. also 
the middle panels of Fig. [Sjl despite the different values of 
Lb- The disk-halo transitions at ii = 4.5 kpc and R = 10.5 
kpc occur at 2 = 1.35 kpc and z = 0.45 kpc, respectively. 
This means that most of the ejecta of the inner fountain 
never reaches the hot halo, while the opposite is true for the 
outer fountain (Fig. [8]). 

Once the expanding phase of the evolution ceases, the 
ejecta of the inner fountain falls back almost entirely to the 
midplane under the action of the surrounding pressure of the 
disk. In the case of the outer fountain, instead, the fraction of 



the ejecta that reached the hot halo falls back rather slowly, 
and "fluctuates" above the disk as in RM. 

As anticipated above, the different dynamical behaviour 
of the ejecta influences its loss of angular momentum. In 
the case of the inner fountain, the ejecta does not interact 
substantially with the hot halo and suffers only a small re- 
duction of its angular momentum; moreover, as the ejecta 
is pushed back to the midplane, it is dragged by the disk 
and its angular momentum increases. Besides the fact that 
the ejecta of the outer fountain emerges from the disk with 
a larger angular momentum (per unit mass), it also loses 
a larger amount of it because its interaction with the hot 
gaseous halo lasts longer and involves a larger fraction of its 
mass. 
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Finally, the bottom panel of Fig. [10] illustrates how the 
amounts of ejecta moving inward and outward (relative to 
the fountain position) change with R. After 150 Myr the 
fraction of ejecta moving inward amounts to 20%, 30% and 
40% at ii = 4.5 kpc, R = 8.5 kpc and R = 10.5 kpc, respec- 
tively (cf. also the thin lines in the middle panel of Fig. 
Thus, this quantity increases with R, in pace with the loss of 
angular momentum. The fraction of the ejecta moving out- 
ward, instead, is lower and is about the same (~ 10%) for 
RM and the outer fountain, but it is comparatively larger 
for the inner fountain, as this one is losing less angular mo- 
mentum. 



5.2 The influence of the gaseous halo 

In the previous sections we have seen how the gaseous halo 
influences the circulation of the fountain gas. In this sec- 
tion we compare the standard model with an analogous one 
where the halo is missing. This does not mean that extrapla- 
nar gas is absent; actually, the density of the ISM decreases 
along the z direction following Eq.[l]and its circular velocity 
depends on z, as illustrated in Fig. [T] 

As shown in the upper and middle panels of Fig. 1111 in 
this model the loss of angular momentum by the ejecta is 
reduced compared to RM (Fig. [5]), although a larger amount 
of the fountain gas is lifted up. In fact, the absence of the 
coronal gas favours the rise of the fountain, but the ejecta 
interacts with an ambient medium which is also rotating (al- 
though at a lower rate than that of the midplane), and the 
transfer of angular momentum is limited. The ambient gas 
density at large z is now higher than in the presence of the 
hot halo. Radiative losses are thus facilitated, and the ejecta 
does not loiter above the disk, as in SMFl, but returns en- 
tirely back to it at later times (cf. the middle panel of Fig. 
Ilip . This explains why, despite a lower loss of angular mo- 
mentum, a larger fraction of ejecta moves inward compared 
to RM (bottom panel). In fact, on the disk the centripetal 
force is larger; as more angular momentum deficient gas falls 
on it, a larger amount of ejecta move inward. Likewise, the 
fraction of ejecta moving outward increases compared to RM 
due to the lower loss of angular momentum that facilitates 
the spreading of the material outward the injection region. 

In conclusion, a smaller amount of gas (40%) remains 
within 7.5 < R < 8.5 kpc compared to RM because a larger 
fraction of it moves both outward and inward. 



Figure 10. Evolution of two models of fountain powered by the 
same number of SNe II as in RM, but located at i? = 4.5 kpc 
(thick lines) and R = 10.5 kpc (thin lines) The meaning of the 
panels, as well as that of the line types is the same as in Fig. |5] 



— 5.3 Neglecting the neutral hydrogen 



It is known t hat the neutral hyd rogen distribution is rather 
clumpy (e.g. IWolfire et al]l2003l l. instead of smooth, as we 
have assumed (cf. Eq. [l]). Actually, such an assumption is 
necessary because an unrealistic grid resolution would be 
needed in order to resolve spatially the H I clouds. One can 
expect that the leading shock wave of the bubble/fountain 
moving through a realistic ISM would advance faster, be- 
cause most of the volume is filled by the rarefied H II gas 
which allows higher shock velocities, as well as lower post- 
shock radiative losses. To evaluate how our assumption on 
the H I distribution influences our results, we have run a 
model similar to the reference one, but where H I is absent; 
this is a rather extreme assumption, but allows us to es- 
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tablish an upper limit value on the heights reached by the 
fountain. 

Actually, after f 50 Myr the amount of the ejecta above 
800 pc and above 1200 pc in the model without H I is almost 
doubled compared to RM (cf. the second panel of Fig. [TT|) . 
This is not a consequence of larger heights reached by the 
ejecta, because it actually attains approximately the same 
altitudes as in RM. In fact, we find that the lifted ejecta does 
not fall back and remains essentially at the height reached 
initially; at larger values of z the amount of ejecta even tends 
to increase with time. This behaviour is due to the decrease 
of the radiative losses as a consequence of the lower density 
of the ambient medium: the hot ejecta does not cool and 
tends to keep moving upward by buoyancy. 

As a larger fraction of ejecta interacts longer with the 
hot halo compared to RM, the loss of angular momentum is 
also larger, as well as the amount of ejecta moving inward (cf. 
the first and third panels of Figs. [S] and [TT] this mechanism 
is substantially similar to that discussed in section [STTI with 
regard to the model with R — 10.4 kpc. 

In conclusion, in the absence of H I, 40% of the ejecta 
mass remains above — 800 pc, and less than 40% is found 
within 7.5 < R < 8.5 kpc (characterizing a larger spread- 
ing)- 

5.4 The role of rotation 



Previous authors (e.g. lde Avillej[2000t ) have investigated nu- 
merically the mechanism of the galactic fountains, focussing 
only on a restricted region of the galactic plane; in their 
computational domain all the quantities (gas and stellar dis- 
tribution, as well as the gravity) depends only on z. In their 
simulations also the Galactic rotation is obviously absent. 
To evaluate the importance of the rotation, we have run a 
model similar to the reference one, but where the rotation is 
inhibited. In order to get a stable disk, we also assumed that 
no radial dependence is present, the ISM being structured 
as a homogeneous disk whose density and pressure depends 
only on z, and the value of these variables at z = 0, as well 
as the vertical scale height, are the same found at J? = 8.4 
kpc in the reference model. The same is true for the ver- 
tical component of the gravitational force, while its radial 
component is set to zero. 

Figure [12] compares results obtained in absence of ro- 
tation with analogous results for the reference model. Both 
models have a maximum resolution of 26 pc. The upper 
panel shows that in the non rotating case the ejecta reaches 
larger heights than in the reference model. In fact, in the 
non rotating model the fountain resembles a column, and 
the SNe II located at its bottom can push the gas upward 
more effectively. In the rotating model, instead, the gas at 
larger z (ejected at former times) lags systematically be- 
hind the freshly expelled gas along the circular trajectory 
(because of the drag by the gaseous halo); the fountain is 
thus less efficient at driving the gas to larger heights. 

The differences in the distribution of the ejecta between 
RM and the non rotating model are illustrated in the lower 
panel of Fig. 1121 Contrary to the previous cases, we now do 
not separate ideally the disk in an inner and an outer region; 
instead, we consider a cylindrical volume extending verti- 
cally and intersecting the disk with a circle with a diameter 
of 1 kpc and centred on the cluster located at _R = 8.4 kpc. 
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Figure 11. Evolution of two models of fountains similar to RM, 
but one without the hot gaseous halo (thick lines), and another 
without H I gas (thin lines). The meaning of the panels as well 
as that of the line types is the same as in Fig. |S] 
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Figure 12. The upper panel shows the amount of the ejecta at 
different heights from the Galactic plane obtained by the RM 
(thin lines) compared with those given by a model without ro- 
tation (thick lines). The lower panel displays the amount of the 
ejecta located inside (dashed lines) and outside (solid linos) the 
cylindrical region above the stellar cluster powering the fountain, 
as well as the fraction located outside the cylinder (see the text 
for the details). 



At early times the ejecta in the non rotating model is mostly 
located inside the cylinder, but then straggles outside while 
falling back. The ejecta in the RM shows the same quali- 
tative behaviour, but the mass fraction inside the cylinder 
remains always higher than in the non rotating case. The 
larger the heights reached by the ejecta, the larger its diffu- 
sion: the fraction left behind in RM along a bow shaped tail 
(cf. Fig. [3]) is negligible and does not play a substantial role 
in this analysis. 



6 DISCUSSION AND CONCLUSIONS 

We studied in detail the development of the Galactic foun- 
tains in order to understand their dynamical evolution, their 
role in building intermediate velocity clouds and condensa- 
tions up to few kpc over the disk, their interaction with 
the hot halo and their influence in the redistribution of the 
freshly delivered metals over the disk. In this paper we de- 
scribe the evolution of a single fountain, i.e. a fountain pow- 
ered by SNe II belonging to a single stellar cluster. In a 
companion paper we will present the results of simulations 
of multiple fountains fuelled by SNe II of different OB asso- 
ciations of various richness, occurring stochastically in space 
and time. 

Not all the SNe II contribute to build galactic foun- 
tains. Only bubbles powered by a luminosity a few times 
larger than the critical luminosity Lb break out of the disk 
pushing the ambient gas some kpc above the Galactic plane 
(cf. I2.2.1|l . As shown in Fig. (9] this critical luminosity de- 
creases exponentially with the Galactocentric distance. In 
the next paper we will show how, combining this behaviour 
with the cumulative probability for the occurrenc e of A^'sn 
clustered supernovae (|Higdon fc Lingenfehej[2005h . it turns 
out that only ^ 60% of all SNe II occurring on the disk can 
give rise to Galactic fountains. This fraction reduces almost 
linearly with radius; for instance, it becomes ^ 40% within 
7? = 4 kpc. It is likely that, in the case of multiple fountains, 
also weaker superbubbles may reach larger heights above the 
plane finding their way through breaks and tunnels carved 
by previous, adjacent and more powerful superbubbles. A 
discussion on this point is postponed in the companion pa- 
per. 

In any case, a non-negligible fraction of superbub- 
bles seems to remain contained by the thickened Galac- 
tic disk and do not participate to a large scale gas circu- 
lat ion. The fate of suc h superbubbles has been envisaged 
bv iTenorio-Tagl j (Il996l ) who suggests that the superbubble 
stops its expansion stalling at an height 2; ~ 1 kpc. The hot, 
low-density, metal-rich gas of its interior waits for the onset 
of radiative cooling. Small overdense inhomogeneities cool 
faster increasing their density and forming droplets which 
start to fall toward the disk. This rain of small (~ 0.1 pc), 
metal-rich cloudets occurs over a kpc-scale volume of the 
galaxy. These cloudets will be successively photo-evaporated 
by the stellar radiation, eventually mixing with the ISM. Al- 
though some objection may be raised about the details of 
this scheme (jRecchi et al.ll200lh . there is no doubt that a 
stalled bubble must deflate as a consequence of the depres- 
surization of its interior due to radiative losses. It is also pos- 
sible that the turbulences of the local medium distort and 
destroy the superbubble, anticipating its disappearance. In 
conclusion, weak bubbles do not affect much the general dis- 
tribution of the ambient medium. Their size can be rather 
small, of the order of 100 pc (cf. the left panel of Fig. ??), 
and their metal content remains strongly localized. 

On the contrary, powerful superbubbles give rise to 
Galactic fountains, and are e xpected to drive a gas circu- 
lation on large scale (IShapiro fc Field 1976. ; .Bregman.l98Q ) . 
Sim ple ballistic simulations de scribed in the Appendix (see 
also lFraternali fc Binnevll2006l ) , where the clouds are treated 
as bullets ejected into the halo from the Galactic midplane, 
actually show that the clouds fall back on the disk at dis- 
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tances that can be as large as several kpc from the starting 
point, depending on the modulus and direction of the initial 
velocity (cf. the Appendix). Our hydrodynamical simula- 
tions, however, show that the majority of the gas lifted up 
by the fountains falls back on the disk remaining within a 
radial distance AT? = 0.5 kpc from the place where the foun- 
tain originated. Our reference model, a fountain powered by 
100 SNe II occurring at _R = 8.5 kpc, may eject material 
up to ~ 2 kpc which than collapses back mostly in form of 
dense, cold clouds and filaments. Although the simple bal- 
listic model indicates that these cloud should move mostly 
outward (cf. the Appendix), our simulations show instead 
that they have the tendency to drift toward the Galactic 
centre as they approach the disk. This effect is due to the 
interaction between the fountain and the gaseous halo above 
the disk. This halo is at rest and exerts a drag force on 
the rotating fountain gas which looses part of its angular 
momentum and is pulled inward by the gravitational force. 
Nearly 60% of the metals delivered by the SNe II remains in 
the radial range 8 < i? < 9 kpc, while the rest is dispersed 
on a larger range, mostly inward rather than outward (cf. 
Fig. [5]). We stress that even in absence of the gaseous halo 
the fraction of metals moving inward still prevails over that 
moving outward, although to a smaller extent (cf. Fig. fTTj) . 
In fact, the exponential tail of the vertical distribution of 
the disk gas is dense enough to exert a non-negligible drag 
on the fountain. 

We point out that ou r results are different from those 
of iBooth fc Theund (120071 ') who obtain much longer trajec- 
tories for the clouds pushed by the fountains. As an exam- 
ple, these authors find that the gas ejected at 7i = 8, 5 kpc 
ne xt passes through the di sk at _R = 25 kpc (cf. their Fig. 
5'l. lBooth fc Theund l|2007l ') adopt an SPH code and clouds 
are modelled using 'sticky particles' that move ballistically 
through the ambient gas. However, also in the framework 
of pure ballistic trajectories (which is the most favourable 
to obtain large ra dial displacements) it is difficult to recon- 
cile the results of iBooth fc Theunj (120071 ) with the shorter 
trajectories illustrated in our Appendix. 

We also stress that our results about the radial displace- 
ment of the metals are in agreement with recent chemical 
models of the Milky Way aimed to explain the chem ical 
gradient along the Galactic disk (|Cescutti et al.l [20071 ). In 
these models the disk is approximated by several indepen- 
dent rings, 2 kpc wide, without exchange of matter between 
them. Because such models are successful in reproducing the 
observed gradients, one can infer that most of the metals do 
not move radially more than 1 kpc from the place where 
they are created, as found by our simulations. 

We finally make some remarks about the interaction be- 
tween the fountains and the hot gas of the halo. The SN II 
feedback, infact, is an important element to understand the 
energy budget of the halo (e.g. Wang 2005) . Unfortunately, 
at this stage, a definite conclusion can not be drawn by our 
models of single fountains. The amount of thermal energy of 
the gas of the halo in our Galactic model is Eth ~ 10^^ erg. 
Although this gas is set initially at rest, some angular mo- 
mentum is transferred to it by the rotating ISM of the disk 
through numerical viscosity. While the amount of this angu- 
lar momentum is small (and the induced rotation is highly 
subsonic), it introduces some noise at the vertical edges of 
our computational box, and the thermal energy of the gas 



reduces of an amount AEth = 3 x 10^^, about 3% of the 
total. This energy must be compared with Etnt = 5 x 10^^ 
erg, which is the energy injected in the numerical grid by 
the SNe II. This energy is ten times lower than the per- 
turbation generated at the boundaries, and any evaluation 
of the energy transfer from the fountain to the halo is pre- 
vented. This subject will be toroughly investigated in the 
companion paper on the multiple fountains, when the en- 
ergy injected by the SNe is ~ 10^* erg and we will be able 
to give a reasonably accurate evaluation of the fraction of 
energy transferred to the halo gas. 
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APPENDIX A: THE BALLISTIC CASE 

To compare the hydrodynamical description of a single foun- 
tain give n in this paper with the pure kinematics picture 
given bv lBreemanl (fl980h . we have built a model to follow 
the orbits of gas clouds that are ejected from the disk and 
move ballistically within the same potential well adopted in 
our hydrodynamical models (see sect. \2.l\ . 

In this simple model we integrate the trajectory of only 
a single cloud which is modelled as a particle initially moving 
around circular orbits in the plane z = and receiving a kick 
velocity aloft. The integration is performed in cartesian co- 
ordinates using the fourth-order Runge-Kutta method, and 
the trajectory is followed until the cloud falls back to the 
disk. 

In all the cases presented in Fig. lAll the cloud starts 
at _R = 8.5 kpc with a kick velocity v^] = lO'^ km s~^. 
The cases differ only in the direction of Vcj. In the bottom 
panel of Fig. lAll three trajectories in the z — R plane are 
shown for which Vej is perpendicular to the circular velocity 
Vc, but have three different inclinations with respect to the 
Galactocentric distance R. The three trajectories in the top 
panel refer to cases where Uoj is orthogonal to R but with 
different slope angles to Hoj . The case of orthogonal ejection 
velocity to the Galactic plane is represented in both panels 
for an easier comparison. 
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Figure Al. Trajectories of a single cloud in the z — R plane, all starting at i? = 8.5 kpc with a kick velocity v^j = 10^ km s ^ . See the 
text for the details. 



Figure IXTI emphasizes the importance of the initial con- 
ditions. In the cases illustrated in the bottom panel, despite 
the very different orbits the cloud falls back on the plane 
essentially at the same Galactcentric distance, not far from 
the starting point. In the examples shown in the top panel, 
instead, the cloud can "land" much closer to, or quite fur- 
ther away, the Galactic centre. In these cases one can expect 
consequences in the kinematics of the thick disk, as well as 
in the stellar abundance gradients along the disk. 



